Scaling of geometric phases close to quantum phase transition in the XY chain 
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We show that geometric phase of the ground state in the XY model obeys scaling behavior in 
the vicinity of a quantum phase transition. In particular we find that geometric phase is non- 
analytical and its derivative with respect to the field strength diverges at the critical magnetic field. 
Furthermore, universality in the critical properties of the geometric phase in a family of models 
is verified. In addition, since quantum phase transition occurs at a level crossing or avoided level 
crossing and these level structures can be captured by Berry curvature, the established relation 
between geometric phase and quantum phase transitions is not a specific property of the XY model, 
but a very general result of many-body systems. 
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The phase factor of a wave function is the source of 
all interference phenomena and one of most fundamen- 
tal concepts in quantum physics. Recent considerable 
interest in this field is motivated by the pioneer work of 
Berry Jj. Berry discovered that a geometric phase, in 
addition to the usual dynamic phase, is accumulated on 
the wave function of a quantum system, provided that 
the Hamiltonian is cyclic and adiabatic. Since then, the 
adiabatic geometric phase and its generalizations 0, |^ 
have found many applications to broad fields 0, 0] , such 
as condensed matter physics 0,013, atomic, molecular, 
and optical physics, and quantum computation^, etc.. 

Very recently, CaroUo and Pachos demonstrated the 
close relation between geom etric phases and quantum 
criticality of spin chains "lill. In particular they showed 
that a nontrivial geometric phase difference between the 
ground state and the first excited state exists in the XX 
model if and only if the closed evolution path circulates 
a region of criticality. Quantum phase transitions usu- 
ally occur for a parameter region where the energy lev- 
els of the ground state and the excited state cross or 
have an avoided crossing, and is certainly one of the ma- 
jor interests in condensed matter physics [TlL IT3l | and 
quantum information^^. Geometric phase, as a mea- 
sure of the curvature of Hilbert space, can reflect the en- 
ergy level structures and can capture certain features of 
quantum phase transitions. However, at least two impor- 
tant problems need to be addressed, (i) The XY model 
is parameterized by 7 and A (see the definitions below 
Eq-O)- Two distinct critical regions appear in param- 
eter space: the segment (7, A) = (0,(0,1]) for the XX 
chain and the critical line Ac = 1 for the whole fam- 
ily of the XY modeim Il5| . The nontrivial geometric 
phase difference between the ground state and the first 
excited states calculated in Ref. ^3 can be used as a 
measure of the presence of the first critical region, but 
this measure is unlikely to remain valid for the second 
critical region. Whether one can reveal the latter crit- 
ical region using geometric phase is of significance, (ii) 
As also noted in Ref.^^l) ^ challenging but also impor- 
tant question warranting further study is whether the 
typical features of quantum criticality, such as the scal- 
ing feature, critical exponents and universality, etc., have 



relation to geometric phases in this many body system. 
Answering these questions is certainly significant for a 
deeper understanding of quantum phase transitions, and 
also from the perspective of geometric phases. So further 
results that bridge these two interesting areas of research 
are of great relevance. 

In this paper, instead of using the geometric phase dif- 
ference between the ground state and the first excited 
state as a signature of quantum criticality, we focus on 
the relation between geometric phase of the ground state 
and quantum criticality in the XY chain. We analyse ge- 
ometric phases near the critical point of the XY model 
and find that the geometric phase is non-analytical and 
its derivative with respect to the field strength A diverges 
at the critical line described by Ac = 1. In particular the 
geometric phase obeys scaling behavior in the vicinity of 
a quantum phase transition. Furthermore, universality 
in the critical properties of geometric phase for a family 
of models is verified. These results show that the key 
ingredients of quantum criticality are present in geomet- 
ric phases of the ground state. In addition, we show that 
the relation between geometric phase and quantum phase 
transitions established here is not model dependent, but 
is valid in a wide variety of systems. 

The system under consideration is a spin- 1/2 XY 
chain, which consists of N spins with nearest neighbor 
interactions and an external magnetic field. The Hamil- 
tonian of the system is given by 
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where M = {N — l)/2 for N odd and cr^ (/i = x, y, x) are 
the Pauli matrices for the jth spin. We assume periodic 
boundary conditions. The parameter A is the intensity of 
the magnetic field applied in the z direction, and 7 mea- 
sures the anisotropy in the in-plane interaction. This XY 
model encompasses two other well-known spin models: it 
turns into transverse Ising chain for 7 = 1 and the XX 
(isotropic XY) chain in a transverse field for 7 = 0. 

As for quantum criticality in the XY model, we need 
to distinguish two universality classes depending on the 
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anisotropy 7. The critical features are characterized in 
term of a critical exponent v defined by ^ ~ |A — Ad""^ 
with ^ representing the correlation length. For any value 
of the anisotropy 7, quantum criticality occurs at a criti- 
cal magnetic field Ac = 1. For the interval < 7 < 1 the 
models belong to the Ising universality class character- 
ized by the critical exponent v = 1^ while for 7 = the 
model belongs to the XX universality class with v — \/2 

To investigate the geometric phase in this system, we 
introduce a new family of Hamiltonians that can be de- 
scribed by applying a rotation of around the z di- 
rection to each spin, i.e., = g^Hg]^ with = 

. J exT>(i6af /21 . The critical behavior is inde- 
pendent of (j) as the spectrum (see below) of the system 
is (f) independent. This class of models can be diagonal- 
ized by means of the Jordan- Wigner transformation that 
maps spins to one-dimentional spinless fermions with cre- 
ation and annihilation operation aj and via the rela- 
tions, aj = ([li^jCrf)aj 0,0. Due to the (quasi) 
translational symmetry of the system we may introduce 
Fourier transforms of the fermionic operator described by 
dk = ^X]jaiexp(-i27rjfc/7V) with k = -M, • • • , M. 
The Hamiltonian can be diagonalized by transform- 
ing the fcrmion operators in momentum space and then 
using the Bogoliubov transformation. The result is given 
by H — X]fe^'£(clcfe — 1), where the energy spectrum 

Ak = y^(A - cos(27rfc/iV)) + 72 sm'^{27rk/N) and Ck = 
dk cos ^ — id^^e^*"^ sin ^ with the angle 9^ defined by 
cos^fe = (cos2^- A)/Afc. 

The ground state \g) of is the vac- 

uum of the fermionic modes described by 
Ck\g) — 0, and can be written as \g) — 
UtLi (cos^|0)fc|0)_fc-ze2^*sinf |l)fc|l)_fc) , where 
|0)fc and |l)fc are the vacuum and single excitation of the 
kth mode, d^, respectively. The ground state is a tensor 
product of states, each lying in the two-dimensional 
Hilbert space spanned by |0)A:|0)_fc and |l)fe|l)_fc. 
The geometric phase of the ground state, accumulated 
by varying the angle from to tt, is described by 
~ ~JT lo i9\^4>\9)^'t'-^^\' ^^'^ found to be 



M , 



fc=i 



cos(27rfc/iV) - A 
Afc 



(2) 



To study the quantum criticality, we are also interested 
in the thermodynamic limit when the spin lattice number 
TV -^00. In this case the summation X]fc=i 
be replaced by the integral i JJ^ dip with if> — ; the 
geometric phase in the thermodynamic limit is given by 
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cos 93 — A 



\l (cos Lp — A)2 -I- 72 sin^ Lp ^ 



dip. (3) 



To demonstrate the relation between geometric phase 
and quantum phase transitions, we plot geometric phase 
(3g and its derivative d(3g/d\ with respect to the field 
strength A as a function of the Hamiltonian parame- 
ters A and 7 in Fig. 1. Two particular features are 
notable: (i) The non-analytic property of the geomet- 
ric phase along the whole critical line Ac = 1 in the XY 
model is clearly shown by anomalies for the derivative 
of geometric phase along the same line; (ii) Geometric 
phase for the XX model under the thermodynamic limit 
is very special in the sense that, instead of using the geo- 
metric phase difference between the ground state and the 
excited phase as the signature of phase transitionjl^], the 
geometric phase of the ground state itself also serves the 
same role. Geometric phase under the thermodynamic 
limit can be obtained explicitly from Eq.Q for 7 = as 



/3g 



27r — 2 arccos(A) when A < 1 and fig = 27r when 



A > 1. Thus, a nontrivial geometric phase of the ground 
state itself is also a witness of quant um p hase transition; 
this fact has also been shown in Ref . . 




FIG. 1: (color online), (a) Geometric phase j3g of the ground 
state (b) and its derivative dj3g/d\ as a function of the Hamil- 
tonian parameters A and 7. The lattice size A'' = 10001. 
There are clear anomalies for the derivative of geometric phase 
along the critical line Ac = 1. 

To further understand the relation between geometric 
phase and quantum criticality, we investigate the scaling 
behavior of geometric phases by the finite size scaling 
approach [1^. We first look at the Ising model. The 
derivatives d(3g/d\ for 7 = 1 and different lattice sizes are 
plotted in Fig[21 There is no real divergence for finite N , 
but the curves exhibit marked anomalies and the height 
of which increases with lattice size. The position An, of 
the peak can be regarded as a pseudo-critical point 
which changes and tends as iV~^-®°'^ towards the critical 
point and clearly approaches Ac as A'^ ^ 00. As shown in 
Fig. 3(a), the value of dj3g/d\ at the point Am diverges 
logarithmically with increasing lattice size as: 



« Kiln at const. 

d\ 



(4) 



with Ki = 0.3121. On the other hand, as shown in Fig. 
3(b), the singular behavior of d(3g/d\ for the infinite Ising 
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FIG. 2: (color online). The derivatives dPg/dX for the Ising 
model (7 = 1) as a function of the Hamiltonian parame- 
ter A. The curves correspond to different lattice sizes A'" = 
21,101,501,1001,00. With increasing the system sizes, the 
maximum becomes more pronounced. The inset shows that 
the position of the maximum changes and tends as 
towards the critical point Ac = 1. 



chain can be analyzed in the vicinity of the quantum 
criticality, and we find the asymptotic behavior as 



K2 In |A — Ac 



const. 



(5) 



with K2 = —0.3123. According to the scaling ansatz in 
the case of logarithmic divergence the ratio |k2/'«i| 
gives the exponent v that governs the divergence of the 
correlation length. Therefore, ~ 1 is obtained in our 
numerical calculation for the Ising chain, in agreement 
with the well-known solution of the Ising model [TU l. 
Furthermore, by proper scaling and taking into account 
the distance of the maximum of pg from the critical 
point, it is possible to make all the data for the value 
of F = [1 — exp((i/3g/dA — dPg/dX\\^)] as a function of 
7V^/%A — Ac) for different N collapse onto a single curve 
[1^ . The result for several typical lattice sizes in the 
Ising model is shown in Fig. 4, where we can also extract 
the critical exponent v = 1. 

A cornerstone of quantum phase transitions is a uni- 
versality principle in which the critical behavior depends 
only on the dimension of the system and the symmetry 
of the order parameter. The XY model for the interval 
7 e (0, 1] belong to the same universality class with crit- 
ical exponent i' — 1. To verify the universality principle 
in this model, we check the scaling behavior for different 
values of the parameter 7. The asymptotic behaviors are 
also described by Eqs. Q and For instance, from 
Fig. 3 we get ki ~ 0.5234 and K2 = -0.5238 for 7 = 0.6. 
Moreover, we also verify that, by proper scaling, all data 
for different N but a specific 7 can collapse onto a single 
curve. The data for 7 = 0.6 are show in Fig. 4. We can 
extract the same critical exponent = I from the data 
shown in both Fig. O and 01 

Comparing with the 7 7^ case, the nature of the di- 
vergence of d(3g/dX at 7 = belongs to a different univer- 
sality class, and the scaling behavior of geometric phase 
can be directly extracted from the explicit expression of 
geometric phase fig = 2-k — 2arccos(A) (A < 1) in the 
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FIG. 3: (color online), (a) The maximum value of the deriva- 
tive dPg/dX at the pseudo-critical point Am as a function of 
lattice sizes. The slope of the line is 0.3121 (0.5234) for 7 = 1 
(7 = 0.6). (b) The derivatives dj3g/d\ for the thermodynamic 
limit logarithmically diverge on approaching the critical value. 
The slope of the line is -0.3123 (-0.5238) for 7 = 1 (7 = 0.6). 
The ratio between the two slopes in (b) and (a) for a fix pa- 
rameter 7 is the critical exponent v. Here ~ 1 is obtained 
for both 7 = 0.6 and 1, as expected by the concept of univer- 
sality in the XY model. 
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FIG. 4; (color online). The value of = [1 - exp(d/3g/dA - 
dPg / d\\\^)] as a function Y(A — Ac) for different lattice sizes 
N = 51, 101,501, 1001. All the data for a fix parameter 7 
collapse on a single curve, as expected from the finite size 
scaling ansatz. 



thermodynamic limit. Since df3g/d\ — \/2{l — A)^^/^ 
(A ^ 1~), we can infer the known result that the criti- 
cal exponent i/ = 1/2 for the XX model. Therefore, the 
above results clearly show that all the key ingredients 
of the quantum criticality are present in the geometric 
phases of the ground state in the XY model. 

We now demonstrate that the relation between geo- 
metric phase and quantum phase transitions addressed 
above is valid in a general case: quantum phase transition 
occurs at level crossings or avoided level crossings, and 
these kinds of level structures usually can be captured 
by the geometric phase of the ground state. Consider a 
generic system described by the Hamiltonian H(ri) with 
r] a dimensionless coupling constant. For any reason- 
able r], all observable properties of the ground state of H 
will vary smoothly as 77 is varied. However, there may 
be special points denoted as rjc, where there is a non- 
analyticity in some property of the ground state at zero 
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temperature, rjc is identified as the position of a quan- 
tum phase transition. Non-analytical behavior gene rally 
occur at level crossings or avoided level crossings On 
the other hand, we also consider geometric phases in a 
generic many-body system where the Hamiltonian can 
be changed by varying the parameters R on which it 
depends. The state \tp{t)) of the system evolves accord- 
ing to Schrodinger equation iJ(R(t))|-0(i)) = ihdt\ip{t)) . 
At any instant, the natural basis consists of the eigen- 
states \n(R)) of iJ(R) for R = R{t), that satisfy 
i7(R)|n(R)) E„{R)\n{R)) with energy E„{R) (n = 
1, 2, 3 • • • ). Berry showed that the geometric phase for a 
specific eigenstate, such as the ground state {\g) = |1)) 
of a many-body system we concern here, adiabatically 
undergoing a closed path in parameter space denoted by 
C, is given by /3g(C) = - / Vg{R) ■ dS, where dS de- 
notes area element in R space and Vg{R) is the Berry 
curvature given byQ 



VgiR)=ImJ2 



{g\\/-EtH\n){n\\/RH\g) 



n#9 



[En 



(6) 



The energy denominators in Eq.© show that the Berry 
curvature usually diverges at the point of parameter 
space where energy levels are cross and may have maxi- 



mum values at avoided level crossings. Thus level cross- 
ings or avoided level crossings, the two specific level 
structures related to quantum phase transitions, are re- 
flected in the geometry of the Hilbcrt space of the system 
and can be captured by geometric phases of the ground 
state. Therefore, the relation between geometric phase 
and quantum phase transitions demonstrated herein is in 
fact a very general result and not a specific property of 
the XY model. 

In summary, we established the close relation between 
geometric phase of the ground state and quantum phase 
transitions in a general many-body system. As a typical 
example, we show in detail that all the key ingredients 
of quantum criticality, such as scaling features, critical 
exponents and universality, etc., are present in the geo- 
metric phases in the XY spin chain. 
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reading of this paper. This work was supported by NSF 
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Note added- After this paper was completed, I got a 
manuscript 16] where the general relation among Berry 
phases, topology, and quantum phase transitions in many 
body systems was independently established. 
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